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Abstract 

We calculate next to leading order QCD corrections to the direct 
decays h — > J/^' + X and b — > r}c + X. The strong renormalization 
scale dependence is seen to persist also in this order and in fact the 
rate is driven to an unphysical negative value at the Mi, scale. We 
show that this is a consequence of the strong suppression and scale 
dependence in the leading order term. Large cancellations take place 
between terms from three different orders in a^. Lacking a third order 
calculation we are forced to throw away the leading order term and all 
but one of the next to leading order terms. The remaining cc gluon 
emission term could very well be the dominant one. Even if this is not 
the case, the picture of a rate at least as suppressed as in the leading 
logarithmic approximation survives. 



Recently, there has been renewed interest in the theoretical treatment 
of various decays involving B mesons (see e.g. ^ ||, H, ^). Since it now 
appears hkely that dedicated B "factories" may be built in the near future, 
it is of great importance to match the new level of experimental accuracy 
to state-of-the-art calculations within the Standard Model. Only then can 
one gain sensitivity to eventual new physics in the study of B decays. Even 
within the Standard Model, one has to make sure that the strong interaction 
physics involved when relating the Kobayashi-Maskawa origin of CP violation 
to measured B decays is well enough understood to establish this final piece 
of the six-quark Standard Model. 

Presently, most of the B decay phenomenology seems to be reasonably 
well understood, with some notable exceptions. First, the measured semilep- 
tonic branching ratio of B mesons seems to be too low (i.e., the non-leptonic 
branching ratio seems to be too high) compared to theoretical calculations. 
This has already led the authors of |Q to put the question whether there is 
new, exotic physics involved in B decays. Secondly, the QCD-improved ef- 
fective non-leptonic Hamiltonian seems to give a too small a branching ratio 
of i? — > J/\l/ + X, even if cascade decays are taken into account 

Both the semileptonic and nonleptonic decay rates have been calculated in 
the next to leading logarithmic approximation. In the calculation of the next 
to leading order correction to the nonleptonic rate, quark masses were put to 
zero. After phase space integration this leads to a complete cancellation of 
one of three next to leading order terms (each with different Wilson coefficient 
structure). It has been argued that neither higher order QCD corrections nor 
nonperturbative corrections are large enough to explain the low semileptonic 
branching ratio [Q. 

In the case of the B J/\E' + X decay, a next to leading order calculation 
of the perturbative QCD corrections has been lacking (in a previous attempt 
[^] the renormalization group was not treated properly [0]). 

In this paper, we present the results of a calculation of -B J/"^ + X 
in the next to leading logarithmic approximation (for details, see [^). In a 
straightforward way we also extend this to rjc production. 

As will be seen, the reason that the rate for direct S — > J/\l/ + X decay 
comes out small is that the Wilson coefficient of the colour singlet effective 
four fermion operator is strongly suppressed by the QCD evolution and be- 
comes zero just below the Mf, scale. In the next to leading order, colour octet 
contributions from cc bremsstrahlung appear. These are not suppressed by 
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the QCD evolution and could therefore be as important as the leading order 
term. This is obscured in the next to leading order calculation by the strong 
scale dependence of the one-loop contributions, making the rate negative at 
the Mfe scale. We show that this is a consequence of the strong suppression 
and scale dependence of the colour singlet Wilson coefficient, and propose a 
slight modification of the resummation to overcome this problem. After a 
careful analysis, the next to leading order calculation is seen to indicate a 
further suppression of the rate as compared to the leading order at the Mb 
scale. This reinforces the suspicion that the bulk of J/\I's produced in B 
decays originate from other mechanisms. 

The nonleptonic effective Hamiltonian describing the B decays we treat 
here can be written 0, |^ 



Here V are Cabibbo-Kobayashi-Maskawa fiavour mixing matrix elements (we 
neglect the Cabibbo suppressed channel c — *• m compared to c — > s), is 
the Fermi weak coupling constant and N is the number of colours ( three for 
QCD). The O's are four fermion shortdistance operators, all with the same 
{V ~ A) ^ {V — A) 7 matrix structure but with different colour structures, 
while the C (/i)'s are the corresponding renormalization scale dependent Wil- 
son coefficient functions. In the last line of Eq. (|lD the effective Hamiltonian 
is written as a linear combination of a colour singlet (Oi) and a colour octet 
(Og) operator, which is the most relevant form for our process. The operators 
Oj- are introduced for convenience since they do not mix under evolution. 
We give explicitly 
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0± = (02±Oi)/2. (2) 
Eqs. (|I|) and @ lead to the following relations for the Wilson coefficients: 



C2{fi) = (C+(/x) + C_(/x))/2. (3) 

In the leading logarithmic approximation, the C± were calculated long 
ago D^, H, with the result 



where the anomalous dimensions d± are given by 

(0) 
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with 
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(0) ^ ±12 (NT I) 
'± 2N ^ ^ 

and 

/?o = y AT - ^riF, (7) 

where np is the number of active flavours at the scale /i. 
This leads to the decay rate for B ^ J/'^ + X ^ 



^B-.j/i,+x - ^0 Go, (8) 

where Lq and L2 are related to L± in the same way as Cq and C2 are related 
to C± (eq. 1^), and with 

Go = lC{l-xf (1 + 2x), (9) 
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where 



^ = M-- 



r2 



(11) 



Here the nonperturbative couphng of the cc pair to the bound state is 
parametrized by the S wave function at the origin, i?5(0), which can be 
measured (including order corrections) in J/\l/ decays to lepton pairs. 

Due to the strong suppression of Lq in the expression for Fq (Eq. (^), 
the lowest order width is very small, and strongly dependent on the scale 
parameter ^ (Fig. 1(a), 2). In fact, with the choice yU ~ 2.3 GeV , it can 
even be made to vanish. The hope has been that the next to leading order 
corrections would both increase the rate and make it less sensitive to /i. In 
next to leading order, several new contributions to the Wilson coefficient 
functions appear, summarized by the expression 

C±M = M (l + ^B,) (l + Hii^MzMOfl,) . (12, 

We expect the contribution from penguin induced operators to be small 
for our process (for a calculation of their size for the general Ai? = 1 Hamil- 
tonian, see [|ll], |l^). Note also that for colour singlet cc production penguin 
graphs do not appear until second order in a^. We work in the 't Hooft- 
Veltman regularization scheme [|1^, where B± and R± are related by 

R± = B± + (-21 ± 57/iV ± 23iV t 4n^) - (rf±//5o)M, (13) 

(see for a more complete discussion, including the effects of altering the 
regularization scheme) The beta function coefficient /5i is given by 

ft = _ ]^Mn, - ^-^^n,. (14) 



We find, in agreement with 



±B{NTl) 

B± (15) 
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where B = 7 in the 't Hooft-Veltman regularization scheme. After a cal- 
culation of one-loop and bremsstrahlung contributions we find the rate for 
6 — >^ J/\E' + X in the next to leading logarithmic approximation, 



2 , as (Mw) - as (/i) 



- -^0 Co H — CpGo X 

x(2LoL,iR^-R_) + Ll(^+ 



\N-1 N+1 
+^Mcp (^L,'G, - LoL, (G2 - 2BGo) - Lo'G^) . (16) 

Here Cp = [N"^ — 1)/2N. stems from real gluon emission, whereas Gi 
and G2 are a mixture of real and virtual contributions. These functions Gi, 
z = 1, 2, 3 are given by 

Gi = /C (2(1 - xf (5 + 4 x) log(l -x) + 

4a; (1 - 2 x) (1 + x) log(x) + (1 - 2 x) {I - x) (3 + 5 x)^ + 

Go (4 7rV3 + 4 log(l - x) log(x) + 8 Ct2{x)) (17) 

^ /2 (1 - x) (34 + 23 X - 51 + 16 x^) 
= ^[ 2^^ + 

32(l-x)^x, 8(l-x)^ (3-x2) , 

—^^ — log 2 - J \, ^ log 1 - X + 

2-x (2-x) 

4x2 (26- 19 x + 4x2) , ^ \ ^ / , \ ^ ^ 

^ ^ log(a;) + Go 12 log(^) - 4 (18) 



2-x ""^■^V ' V "'^b' 



G3 = /C(^-(l-x) (l + 37x-8x2)--(l-6x)log(x)), (19) 
where Ci2{x) is the dilogarithmic function 

a^ix) = - r dt'-^^^^. (20) 

JO t 



An important term to notice in (|18|) is the scale dependent term proportional 
to 12Go(x) log(/iVM2) in G2. 

The main difference between our result in Eq. ([T6|) and that of |^ lies in 
the treatment of the large logarithmic corrections that come from the box 
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diagrams, and the way they are summed using the renormahzation group. In 
[^] the J/\l/ particle was treated as being a fundamental, colour singlet field 
produced by the leading logarithmic effective Hamiltonian, to which gluonic 
corrections were applied. This does not seem to be correct, since all short 
distance interactions that may eventually give rise to a J/\l/ should be added 
coherently. There could possibly be other production mechanisms, e.g., a 
{v/cY suppressed (but not Lq suppressed) soft gluon induced fragmentation 
process. This does not, however, effect the calculation of the hard subpro- 
cess rates. Only after the formation time ~ can the J/\l/ particle be 

considered to be a fundamental, colour singlet state. Indeed, as pointed out 
in [H], the renormahzation group summation of the leading logarithmic terms 
in does not reproduce the 0{a1) perturbative result correctly. Already at 
the outset of the calculation in cancellations are performed between terms 
that in a correct treatment should be multiplied by different combinations 
of Wilson coefficient functions. It is therefore impossible to reconstruct the 
correct result from this calculation. Going the opposite way, trying to re- 
produce numerical values given in for some classes of diagrams from our 
calculations, we still find some minor discrepancies. 

In Fig. 1 (curves (a) and (b)) we show the leading and our next to leading 
order result for S — > J/\l/ + X. As can be seen, the next to leading order 
corrections are very large and the result still depends dramatically on the 
scale /i. For large values of /i the rate is even driven to unphysical negative 
values, suggesting that higher order corrections are still anomalously large. 

To understand this behaviour we look at the expansion of the leading 
order term GqLqI^"^) around some scale /i* close to the zero point of Lq 



L^(/i)Go = Lg(/i*)Go + ^^^^^12GoLo(/i*)L2(/i*)log(/iV/x*^) + 

(^^fiin^zj {SGGoLli^n log^(/xV/i*^) + O (Lo(/i*)))(21) 

We note that the renormahzation scale dependence is dominated by the sec- 
ond order term in this expansion. In the next to leading order (|T^), the first 
order term in (pT]) is cancelled by the /i dependent term in G2 (p^: 
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-^iM^12GoLo(/x)L,(/.)log(/.V/i*^) = 
«s(/i*)C^_^2GoLo(/.*)L2(/i*) log(/iV/i*') - 
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(^'^S^)V2Goi^^(^*)log^(/iV-"*') + OiLoif.*))) (22) 

After the cancellation only a second order fi dependent term remains. This 
term is twice as large as the one dominating the renormalization scale depen- 
dence of the leading order term (^I]) , however, and moreover has the opposite 
sign. 

To conclude, the Lq^ suppression postpones the cancellation of the leading 
order renormalization scale dependence until the second order. The next to 
leading order corrections merely mirror the leading order scale dependence. 
Thus, if the appropriate scale for the process is close to the zero point of 
Lq (curve (a) in Fig. 1) the qualitative picture of the scale dependence in 
the next to leading logarithmic approximation (curve (b) in Fig. 1) is an 
inevitable consequence of the strong scale dependence and suppression in 
the leading logarithmic approximation. 

To overcome this problem we propose a simultaneous expansion in Lq 
and as- This amounts to a rather modest resummation, mixing terms from 
three orders in as- In this way terms proportional to Lq or Lq^ are added 
together with the terms that cancel their extreme scale dependence. Since 
we have already seen that we do not gain any information by adding these 
terms separately we do not lose anything by doing the expansion this way. 
Let us now analyze if we can gain something. 

In the new L^-as expansion we find that the old leading term Lq^Gq is 
replaced by the square of the two bremsstrahlung diagrams from the cc pair 
— the G3 term in Eq. (0). We note that the Gs term dominates Lq^Gq 
only at scales close to the zero point of Lq. From the argument above we 
know, however, that the extreme scale dependence of Lq^Gq will be cancelled 
in the new scheme, and it is only after this cancellation that we can expect to 
get something numerically small, as compared to the G3 term. The correct 
question to ask is obviously - will the next to leading order term in the Lo-a^ 
expansion be numerically small? 
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In a complete next to leading order calculation in the new scheme one is 
forced to calculate all terms that are not Lq suppressed. Such a calcula- 
tion could possibly be performed since all corrections coming from two-loop 
diagrams or from higher order matching and evolution of the effective weak 
theory are Lq suppressed. 

We have already calculated several of the terms needed in such a new 
next to leading order calculation. Lacking a complete calculation we try to 
get as much information as possible from the calculated terms and from the 
structure of the unknown terms. The missing terms are all proportional to 
a.^L2 and are therefore suppressed as compared to the leading order 
term. We know, however, that they contain the term SGGq log^(Ai^) that 
cancels the scale dependences in the Go and G2 terms (^,^). It is the 
large numerical prefactor of this log^ term that threatens the viability of the 
expansion. 

The G2 term in eq. (0) is dominated by the cross term between a group 
of one-loop diagrams and the born graph. It is the square of this group of 
diagrams that contain the dangerous log^ term. Since the UV-divergences of 
the one loop diagrams have the same gamma matrix structure as the born 
diagram, we can use G2 / (4Go) as a rough estimate of the square of the 
one-loop diagrams. If we add this approximate term to the next to leading 
order terms that we have calculated exactly, we get 



- aT + livTT + iV^J ) + 

Go (^Lq'^ + -^-^CfLqL22Ri + C%L2^Ri^ — 

— ^^Ci? {G2 — 2BGq) (^LqL2 + -^^CfL2^Ri^ + 

where Ri = R+ — R- and Aus = {as{Mw) — ctsifJ^)). 

Numerically we find that the cancellation is quite effective (Fig. 1 (c)). 
Remember that the born term, which is several times larger than the new 
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leading order term at the Mb scale, is one of the terms in the correction. 
This incomplete second order calculation should not be viewed as an approx- 
imation to the complete one. It is only meant to show that the cancella- 
tion between the terms which are the largest, viewed separately, does take 
place. We have completely neglected the double bremsstrahlung and one- loop 
bremsstrahlung diagrams appearing at the as^L2^ level. The next to leading 
order correction in the L^-ag expansion could very well be substantial, but 
there is reason to believe that it is at least not larger than the leading term. 
In Fig. 2 we display the scale dependence of the relevant ingredients in our 
scheme Lq, L2 and ag separately. 

In the numerical calculations we used Aq^;^ = 190 MeV (as(M^) = 
0.115), Mfe = 5 GeV and x = AM^/M^ = (M^/Mbo)^ ^ 0.38. 

Varying Aqcd and (with M^/Mh fixed) we find - not surprisingly - 
that the Lq-Os expansion works best for large ag and low Mb when the zero 
point of Lq is close to the Mb scale. With ^^(M^) = 0.120 the cancellation is 
almost perfect while it is already beginning to break down at as{Mz) = 0.110. 
Furthermore, we note that the hard gluon emission term G3 is more sensitive 
to variations in Mc/Mb or Ms/ Mb (see for expressions to first order in 
Mg/Mb ) than the other terms. Thus a 5% lower value for Mc/Mb results 
in a 20% increase in the leading order term in the Lo-as expansion. If we 
increase Ms/ Mb from to 0.1 corresponding to Mg ~ 0.5 GeV we find a 5% 
suppression of Go while G3 is suppressed by substantial 17%. 

We have calculated next to leading order corrections also for rjc production 
in B-meson decays, and find an almost identical pattern as in the case of J / ijj 
production. To check the scheme dependence of our calculations, we have also 
utilized the naive dimensional regularization (NDR) scheme, with negligible 
differences in the final results. 

To be able to compare our results with experiments we first have to deal 
with the strong overall dependence on heavy quark masses, the Kobayashi- 
Maskawa matrix element Vcb and on Rs{0), the value of the radial wave 
function for the cc pair at zero distance. Using the measured rate (5.36±0.29 
keV) for the electromagnetic decay of the J/ip |]T3| and the next to leading 
order result 




(24) 



we find |i?5(0)|^ / M^^ ^ 0.5 GeV. This must, however, be considered a very 
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rough estimate since the QCD corrections in (^) are extremely large. In the 
leading order we would get |i?5(0)|^ / ^ 0.2 GeV. Hopefully lattice calcu- 
lations will soon give more precise estimates. To cancel the Vet dependence 
and to reduce the very strong Mf, dependence we normalize our results to the 
semileptonic branching ratio. To next to leading order in «<, the semileptonic 
decay rate is 

r,. = %^M,'g{MjM,) (l - ^/(M./M,)) (25) 

where 

g{x) = l-8x^ - 2Ax^ log(x) + 8a;^ - (26) 

and f{x) can be found in tabulated form in |T^. Using Mc/Mb = Mdo/Mbo ~ 
0.35 in (pSD with as at 5 GeV and the measured semileptonic branching ratio 
(10.7 ±0.5)% [0] we find 

" r,,, 1^ r,, ) Go 

~ (0-2%) X 1^ 0.5Gey )[ 0:05 

The measured inclusive branching fraction for ip production in B decays is 
(1.12±0.16)% . Using measured branching fractions for B decays into ip', Xci 
and for the cascade decays ip' ^ ip + X , ip' ^ xa + 7 and Xci ^ + 1, one 
finds the branching fraction (0.71 ± 0.20) [|ll for direct J/ip production in B 
decays, assuming that no other cascade decays give significant contributions 
to the J/ip production. Using the large value 0.5 GeV for |_R5(0)|^ / and 
taking the size of the QCD-suppression to be ~ 0.05 from Fig. 1, we are still 
two and a half standard deviations from the experimental result. 

We thus conclude from our analysis that unless the experimental value 
goes down as more data accumulate, it seems that the decay B ^ J/"^ + X 
is a process that cannot presently be well described by a standard applica- 
tion of perturbative QCD to zeroth order in the relative velocity of the cc 
quarks. We have pointed out that the strong QCD suppression of the colour 
singlet operator give rise to strong cancellations between terms appearing in 
different orders in as- Although we have suggested a method to ameliorate 
this problem, there is room for other QCD effects, subleading in most other 
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cases, that could bridge the gap between the predicted and observed values 
for the branching ratio. It is thus not at all necessary at this point to make 
the conclusion that exotic mechanisms (i.e. beyond the Standard Model) are 
needed to explain the large experimental rate for this type of decay. 

We thank S. Aminneborg and R. Robinett for useful discussions. This 
work was supported by the Swedish Natural Science Research Council. 
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Figure Captions 



Fig. 1. Renormalization scale dependence of T b-^j/^i+x normalized to the naive 
parton model result Go, in various approximations: 

(a) Leading logarithmic approximation (for consistency, the leading 
order result for is used in the definition of L±). 

(b) Conventional next to leading logarithmic approximation. 

(c) Leading order result in the proposed Lq-Os expansion. 

(d) Incomplete next to leading order result in the Lq-Os expansion, as 
given by Eq. (^31). 

Fig. 2. Renormalization scale dependence of the colour singlet and octet Wil- 
son coefficient functions Lq, L2 and of as- 
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